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Abstract 

The a-decay rates for the nuclides i68,i70,i7i,i72,i73,i74,i76 Y b an d 
148,150,152,154q^ have been estimated from transmission probabilities in a sys- 
tematic a-nucleus potential and from an improved fit to a-decay rates in the 
rare-earth mass region. Whereas a-decay of 152 Gd in natural gadolinium is 
a severe obstacle for the use of gadolinium as a low-energy solar-neutrino 
detector, we show that a-decay does not contribute significantly to the back- 
ground in a ytterbium detector. An extremely long a-decay lifetime of 168 Yb 
is obtained from calculation, which may be close to the sensitivity limit in a 

low-background solar neutrino detector. 

PACS numbers: 23.60.+e,26.65.+t,27.60+j, 27.70. +q 
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The solar neutrino problem is an important current subject to be studied in relation to 
the fundamental physics of neutrino oscillations. After the successful measurements of the 
solar neutrinos by Super-Kamiokande 0J, SAGE 0, GALLEX [|], and SNO |J, the problem 
appears to be solved best in terms of oscillation of v e into other neutrino flavors. To arrive at 
the final solution of the solar-i/ problem, a real-time measurement of the z/ e -energy spectrum 
including pp-, 7 Be-neutrinos is now of central interest ||. 

The solar neutrino detection requires an extremely low-background measurement with 
new technologies. Raghavan |J has suggested a flavor-specific detection scheme with low 
thresholds for the real-time detection of solar neutrinos. The neutrino captures z/ e + 176 Yb — > 
e~ + 176 Lu* and u e + 160 Gd — > e~ + 160 Tb* are based on charged current mediated Gamow- 
Teller transitions. For the identification of a neutrino capture event a delayed coincidence 
between a prompt e~ and a 7-ray from the decay of a nanosecond isomer in 176 Lu or 160 Tb is 
used. For both suggested nuclei, 176 Yb and 160 Gd, one finds sufficient Gamow- Teller strength 
at low energies and a nanosecond isomer which is populated in the neutrino capture reaction. 
Typical detection rates of solar neutrinos are of the order of ten counts per year and per ton 
of detector material 0. 

A new solar neutrino detector using Yb materials has been proposed in LGNS at Gran 
Sasso |5|-|§ , and sources of background for such a detector are presently being studied || . 
The a-decay from the material in a solar-i/ detector is a possible serious problem to pollute 
the spectrum consisting of rare neutrino events. For example, it has turned out that Gd 
material is not suitable for a solar- v detector; the 152 Gd isotope decays by a emission with 
a half-life of Ti/ 2 = 1.08 xlO 14 years, and is followed by the subsequent 148 Sm a-decay with 
T1/2 = 7xl0 15 years. The number of 152 Gd atoms included in the natural Gd metal with 20 
tons amounts to 1.5 xlO 26 atoms. The decay rate calculated by the well-known radioactive 
decay formula 

amounts to about 30,000 counts/sec for the 20 ton Gd material. This is a huge decay rate, 



which decreases the reliability of the coincidence events due to a large number of background 
a-signals and due to accidental coincidences even in the case of a narrow time gate less than 
10~ 8 seconds for coincidence. 

Recently, ytterbium has been suggested as material for a solar-z/ detector. The Q-values 
of a-particle emission for all the Yb isotopes are positive, but relatively small. Hence, each 
Yb isotope must have an extremely long half-life although there has been no report on the 
a-decay lifetime in literature. Such half-lifes may be measurable with a large mass detector 
like a solar-z/ detector. If the half-life of some Yb isotope is relatively short, the a-decay may 
cause a significant background in the analysis of the obtained neutrino data. Therefore, it is 
important to estimate the a-decay rates of the Yb isotopes and to compare them with the 
solar neutrino rates. Recently, Gamow- Teller strengths in the inverse (3 transition 176 Yb — > 
176 Lu were studied, and it has been concluded that Yb-based detectors are well-suited for 
the real-time spectroscopy of the sub-MeV solar neutrinos [fFyiQ|. 

There are several parameterizations for a-decay rates available mainly for nuclei in the 
region of Z > 84 JTT|-|r5|]. The predictions from these parameterizations usually show small 
deviations for experimentally available decay rates. However, there often happen large 
discrepancies, up to ten orders of magnitude, between the predictions for extremely small 
decay rates (corresponding to small Q a values). In this report, we wish to discuss the 
a-decay rates of Yb and Gd isotopes in the rare earth nuclei. 



The lifetime of a-decay was first formulated by Geiger and Nuttall JTJ] . The probability 
of a-particle emission, W, can be expressed as 

W= Pa uT, (2) 

where p a is the probability of finding an a-particle in a nucleus, v the frequency of an 
a-particle appearing at the wall of a nucleus, and T is the transmission coefficient for the 
tunnelling effect. The half-life T\/2 can be parameterized [[HJ as 

^ , _ I ry ,1/3 , _ Z <i 



log l0 T 1/2 = a x + log w -^ + a 2 \jZ d A l l 3 + a 3^= , (3) 



with the half- life Ti/ 2 in years, the Q- value for a-decay Q a in MeV, and the mass (proton, 
neutron) numbers A (Z, N) of the parent nucleus and A d = A — 4 (Z d = Z — 2, N d = N — 2) 
of the daughter nucleus. 

Since the data fittings in most of Refs. [[TT| [3~3j] are focused on the a-decay for the nu- 
clei with Z > 80, we tried to fit the a-decay half-lives for nuclei in the rare earth region 
separately The data for the half-lives of the a-decay nuclei in the rare earth region were 
fitted using Eq. (H) with ai, a 2 , and 03 as three parameters. All Q a values have been taken 
from the mass table of Audi and Wapstra The best fitting parameters, a±, a 2 , and a% 
are a± = —34.577, a 2 = —1.2685, and = 1.7336, respectively. The a-decay half-lives are 
found to be well fitted as can be seen from Fig. [L| The x 2 P er degree of freedom has been 
significantly improved. In our 3-parameter fit \ 2 is reduced by a factor of about 4 compared 



to the 2-parameter fit of Ref. [11], and by a factor of about 2 compared to the 2-parameter 



fit of Refs. [fT2| , p~3|1 which is based on the model of Bethe JT^]. These 2-parameter fits are 
given by |11 



log w T 1/2 = (9MZ° d s /jQ a ) - 58.87 



(4) 



and 113,1131 



log w T 1/2 = 1.598(Z d /y/Q a - Z 2 / 3 ) - 27 '.44 (5) 

Tables | and || show the expected half-lives for the Yb and Gd isotopes. The experimental 
half-life of 1.1 xlO 14 years for 152 Gd agrees with the result of the present fitting within a 



factor of 1.4. However, the predictions from the two parameter fittings by Brown JTTJ and 
Refs. [P2],|l3| differ by three orders of magnitude. A similar agreement is also found for 148 Gd 



and Gd. Again there are strong deviations between the predictions mainly given from 



the fits of a-decay rates for nuclei with Z > 84 |Tl|-p^] . Here, one can address the question 
that the descrepancy may come from the deformation effect on a-decay for the rare earth 
nuclei. But, it should be noted that the most nuclei with Z > 84 are deformed as well, and 
the influence of nuclear deformation deos not change the half-life by an order of magnitude 
(see discussion below). 
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Our prediction for the different quasi-stable Yb isotopes varies from Ti/ 2 = 4.5x10 
years for 168 Yb to Ti/ 2 = 6.2xl0 95 years for 176 Yb. The parameterizations from literatures 
are between one order of magnitude ( 168 Yb) smaller |T2| , p!3| (larger ]TT|) to five orders of mag- 
nitude ( 176 Yb) smaller [fj3,|13|] (larger |L1|). Note that no useful prediction for the relevant 
a-decay half-life of 168 Yb can be derived from Refs. |TT|^T3| because their predictions differ by 
more than three orders of magnitude. Of course, the calculated half-lives of 172 > 173 > 174 > 176 Yb 
exceed any measurable range. But, the excellent agreement between the new fit and the 
half-lives from the folding potential calculations (see below) gives us some confidence that 
our predictions are reliable in a very broad range of half-lives (where other extrapolations 
show huge deviations). 

Because of these large discrepancies between the extrapolations of the various parameter- 



izations, we applied the semi-classical model of Gurvitz and Kalbermann |]T8| together with 



the systematic folding potentials of Atzrott et al. [19]. The required nuclear densities were 



derived from measured charge density distributions |2(J. When no measured charge distri- 
bution was available, the measured distribution of a neighboring isotope was used with an 
adjusted radius parameter R ~ A 1 / 3 . This model was already successfully used to calculate 
the a-decay properties of neutron-deficient so-called p-nuclei in 



(6) 



In the semi-classical model, the a-decay width r a is given as yj 



h 2 

PF — exp 
An 



-2 



k(r)dr 



with the preformation factor P, the normalization factor F 

rr 2 dr 



i k(r) 



(7) 



and the wave number k(r) 



«r) = \l^\Qa-V(r) 



(8) 



where /i, Q a , and r« are the reduced mass, the a-decay energy, and the classical turning 
points, respectively. For + — > + s-wave decay, the inner turning point is at T\ = 0. r 2 



varies from about 7 to 9 fm, and r% varies up to about 350 fm for the lowest Q a value. The 
decay width T a is related to the half-life by the well-known relation r a = h\n.2/Ti/2- For 
the odd Yb isotopes, an additional centrifugal potential was used. The preformation factors 
P were estimated from the systematic trend of P shown in Fig. 2 of Ref. |H] leading to 
P = (4.0 ± 0.8) % for the Yb isotopes and P = (10.4 ± 2.5) % for the Gd isotopes. Larger 
P values can be expected around the shell closures N = 82 and N = 126, whereas smaller 
values for P are found in between. The results of this model are listed in Tables | and |J. 
They are in good agreement with our new fits in the rare earth region. Therefore, the newly 
determined a-decay formula can be used for predictions in the rare earth region with very 
limited uncertainties. 



A systematic study of a-decay half-lives p2| using the same model |T8[ but specially 
shaped potentials is in rough agreement with our results for the Gd isotopes (8 th line in 
Table [H]). Unfortunately, the quasi-stable Yb isotopes are not considered in Ref. f2"2fl . 

There are several sources of uncertainties for the presented predictions for the half-lives 
of the Yb and Gd isotopes. Because of the roughly exponential dependence of the tunneling 
probability on the a decay energy Q a , a factor f(AQ a ) has been calculated by which the 
half-lives are shorted (enlarged) if Q a is increased (decreased) by its uncertainty AQ a . This 
factor f(AQ Q ) is given in Tables | and|T|. Typical uncertainties remain within 40 % with the 
exception of 176 Yb where the very low Q a and its 4keV uncertainty lead to an uncertainty 
of f(AQ a ) = 3.61. This uncertainty is common to all predictions in this work and also in 
Refs. |T3-0. 

The uncertainties of the predictions from the 3-parameter fit in this work are smaller than 
in the previous 2-parameter fits because of the improved \ 2 value. The average deviation 
for logioTi/2 in our fit is about 0.5 corresponding to an average deviation of a factor of 3. 
Therefore, a typical uncertainty of a factor of 3 can be expected for our predicted half-lives. 
The uncertainties in the coefficients a%, ci2, and 03 lead to minor unertainties compared to 
the average deviation of a factor of 3. 

In the case of the folding potential calculations, the predicted values depend sensitively 



on the average preformation factors derived from the systematics in Ref. |2T|. Typical 
uncertainties of P are of the order of 20 — 30 % which lead to identical uncertainties for 
the predicted half-lives. Together with the uncertainties of Q a , this leads to an overall 
uncertainty of better than a factor of two for the folding potential predictions. 

It is well-known that nuclei in the rare earth region show strong deformations with 
typical deformation parameters in the order of ~ 0.3 |23| which may lead to additional 
uncertainties. The influence of deformation on the decay probabilities has been analyzed 
recently in Refs. [p4-27|. In Ref. [24|, it is stated that various calculations of the half- 
lives agree within a factor of 4 with the spherical result. In Ref. ||25|| , it is suggested that 
the deformation effect enhances the decay probability by roughly a factor of 2, whereas in 



Ref. the decay probability for deformed nuclei is reduced by a factor less than 2. In 
Ref. f27fl , the deformation effect leads only to small variations of the half-lives. The good 
quality of our fit to the decay data in the rare earth region confirms that the influence of 
deformation on the a-decay half-lives remains relatively small compared to the uncertainties 
that arise from the predictions of our and of previous fits [|TT|-|i"3]]. 

In summary, we estimated the half-lives of a-decay for the Gd and Yb isotopes. Our 
predictions from the new fit in the rare-earth region and from the folding potential calcu- 
lations agree within their uncertainties discussed above, whereas previous fits show large 
deviations especially for low Q a values and long half-lives. The "short" half-life of 152 Gd 
of Ti/2 = l.lxlO 14 years leads to serious background problems in a neutrino detector. The 
estimated half-life for 168 Yb is between 4.5 xlO 24 years and l.lxlO 25 years. If this predicted 
half-life is correct, the LENS detector with a mass of 20 tons of natural Yb will detect 
about 6~14 counts/year for a-particle emission with Q Q =1.951 MeV from 168 Yb decay, be- 
ing negligble compared with the solar-z/ events. This low rate causes definitely no serious 
background problem. However, the longest a-decay radioactivity so far measured in the past 
may be found in a low-background detector designed for a real-time solar-z/ spectroscopy. 
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FIG. 1. Ratio R = Tffi/Ty$ between the calculated half-lives Tffi with Eq. (|) to the 
experimental half-lives T^,^ for nuclei in the rare earth region. The dotted lines connect isotopes 
from the same element. The error bars are only experimental. Additional errors from the uncer- 
tainty of Q a are not included (see last lines of Tables | and [IJ) • The analyzed range of Q-values 
and half-lives is broad: Q-values cover the range from 1.9 MeV ( 144 Nd) to 7.0 MeV ( 175 Hg), and 
half-lives Ti/ 2 cover 26 orders of magnitude from milliseconds ( 158 W, 169 Pt) to more than 10 15 
years ( 144 Nd, 148 Sm, 174 Hf, 186 0s). 
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TABLES 

TABLE I. Isotopic abundance (%) of Yb, a-decay energy (Q a ), the expected half-lives (7i/ 2 ) 
in unit of year, and the uncertainty factor f(AQ a ) from the Q value error. The parameters of 
Eq. (H) were adjusted to the half-lives of the nuclides shown in Fig. [l|. 





168 Yb 


170 Yb 


171 Yb 


172 Yb 


173 Yb 


174 Yb 


176 Yb 


(%) 


0.13 


3.05 


14.3 


21.9 


16.1 


31.8 


12.7 


Qa (MeV) 


1.95075 


1.73764 


1.55895 


1.31030 


0.94586 


0.74009 


0.57102 


AQ a (MeV) 


0.00412 


0.00146 


0.00139 


0.00139 


0.00145 


0.00160 


0.00417 


T1/2 (exp.) 
















T 1/2 (est.) a 


1.1-10 27 


1.4-10 32 


1.6-10 37 


8.5-10 45 


3.0- 10 64 


3.8-10 80 


7.7-10" 


T 1/2 (est.) b 


5.5-10 23 


2.4-10 28 


9.3-10 32 


7.4-10 40 


4.7-10 57 


1.8-10 72 


5.5-10 89 


T 1/2 (est.) c 


4.5-10 24 


4.L10 29 


3.5-10 34 


1.L10 43 


1.5-10 61 


8.4-10 76 


6.2-10 95 


Ti/ 2 (est.) d 


1.L10 25 


6.9-10 29 


4.7-10 35 


1.0-10 43 


2.1-10 61 


1.9-10 76 


6.9- 10 94 


/(AQ Q ) 


1.24 


1.09 


1.10 


1.13 


1.23 


1.40 


3.61 



I) 



a Parametrization of B. A. Brown fry] ]. 

'Parametrization of R. Taagepera and M. Nurmia [12], K. K. Keller and H. Z. Munzel [13]. 
c this work, from fit to rare earth region. 
d this work, from folding potential, using P = (4.0 ± 0.8) %. 
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TABLE II. Isotopic abundance (%) of Gd, a-decay energy (Q a ), the expected half-lives (Ti/ 2 ) 
in unit of year, and the uncertainty factor f(AQ a ) from the Q value error. The parameters of 
Eq. (H) were adjusted to the half-lives of the nuclides shown in Fig. [I] including 148 > 150 > 152 Qd. 





148 Q d 


150 Gd 


152 Gd 


154 Gd 


(%) 


_ 


_ 


0.20 


2.18 


Qa (MeV) 


3.27121 


2.80887 


2.20458 


0.91991 


AQ a (MeV) 


0.00003 


0.00630 


0.00142 


0.00115 


T1/2 (exp.) 


7.5-10 1 


1.8-10 6 


1.1-10 14 


- 


T 1/2 (est.) a 


7.6-10 3 


7.1-10 8 


3.7-10 17 


2.9-10 59 


T 1/2 (est.) b 


2.0-10 2 


4.4-10 6 


1.8-10 14 


6.7- 10 50 


T 1/2 (est.) c 


2.6-10 1 


1.1-10 6 


1.5-10 14 


4.1-10 53 


T 1/2 (est.) d 


6.0-10 1 


1.7-10 6 


1.4-10 14 


4.1-10 52 


T 1/2 (est.) e 


2.6-10 1 


9.0-10 5 


6.7-10 13 




/(AQ a ) 


1.00 


1.18 


1.06 


1.17 



a Parameterization of B. A. Brown |11]. 

b Parameterization of R. Taagepera and M. Nurmia [12|, K. K. Keller and H. Z. Munzel [13]. 

c this work, from fit to rare earth region. 

d this work, from folding potential, using P = (10.4 ± 2.5) %. 



e from Ref. [22 
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